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Abstract. The separation of contributions from multi- and single-meson 
production in the missing mass spectrum of the quasi-free pn — > pnX reaction 
constitutes a challenging task when the reaction is studied close to threshold. This 
is especially true if the resolution of the mass determination is comparable with 
the excess energy and if the investigated signal appears close to the kinematical 
limit. In this article we outline a method which permits the extraction of the signal 
originating from the creation of a single meson without the necessity of conducting 
model-dependent simulations. For the pd — > pnXpspectator reactions, the method 
allows one to combine events corresponding to multi-meson production at various 
excess energies with respect to the pn — > pn meson process, and hence leads 
to an increase of the statistics needed for the determination of the shape of the 
multi-meson background. 

As an example of the application of the method, we demonstrate that the 
evaluation of the data from the pd — > pnXpsp process according to the described 
technique enables one to extract a signal of the pn — > pnrj reaction whose shape 
is consistent with expectations, supporting the correctness and usefulness of the 
method introduced. 
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1. Introduction 

In recent years, after many very successful studies of the close-to-threshold production 
of mesons in proton-proton collisions [1-5] , several experimental groups have extended 
the range of their investigations to proton- neutron scattering [6-16]. The results 
obtained can then be compared with theoretical models of the production of mesons in 
different isospin channels and this reduces significantly the ambiguities of such models. 
However, the experiments and the evaluation of the data from the proton-neutron 
interaction are much more difficult than measurements of proton-proton reactions. 
This is mainly due to having to use nuclear targets as a source of neutrons. In 
order to define the full kinematics of the pn pnX events one needs to measure 
at least one more particle than for the analogous meson production in proton-proton 
scattering. This is because the determination of the four-momentum vector of the 
colliding quasi-free neutron requires the measurement of the non-interacting residue 
of the target nucleus (spectator) or the decay products of the created meson. In 
comparison to the pp — > pp meson reactions, the efficiency for registering the quasi 
free pn pn meson process is reduced by a large factor and, additionally, the 
resolution in the four-momentum of the investigated meson also worsens significantly. 
Therefore, in the evaluation of the data one faces problems of low statistics, which 
are especially important when the analysis demands subtraction of histograms and 
when the expected signal is close to the kinematical limit. Such is the situation when 
the quasi-free pn pn meson reaction is investigated by registering all the outgoing 
nucleons from the pd — > pappnX reaction, using the missing mass technique to identify 
the desired meson events. The subscript sp denotes the spectator proton which, it is 
assumed, does not take part in the reaction and X stands for the single meson or 
multi-meson system. 

An identification of the pn pn meson reaction on an event-by-event basis is 
impossible if only the nucleons are measured. One can, however, extract the number 
of registered pn — !■ pn meson events from the missing mass distribution provided that 
the contribution of the continuous spectrum originating from multi-pion production 
can be subtracted from the signal. 

In this article, we describe a method for the subtraction of a physical multi- 
pion background from the missing mass spectrum for quasi-free pn pnX reactions, 
where the decay products of the produced meson are not detected. Such a technique 
is clearly required in, e.g., the derivation of the correct number of pn pn/q 
and pn pnrj' events from the missing mass spectra obtained by the COSY-11 
collaboration [14,15,17]. 

In the next section we introduce the method in the idealized case where the 
number of events is infinite. This simplification will allow us to focus first on the 
problem of the proximity of the signal to the kinematical limit. After explaining 
the main ideas of the background subtraction, we will generalize the discussion to 
the situation encountered in real experimental conditions. Thus, in section 3 we will 
demonstrate how in the case of the low statistics one can use all events taken at 
negative excess energies and combine them to construct the shape of the background. 
Section 4 will be devoted to the justification of the assumptions made and finally in 
section 5 we will show an example of the application of the technique to data taken 
at the COSY-11 facihty. 

Though the method was developed for the analysis of the COSY-11 pd pnijpsp 
results, the conclusions are valid also for the production of the t]' and other narrow 
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mesons. Moreover, since the technique is independent of the detection system, it may 
be also appUed when evaluating data obtained at other facilities. 

2. Derivation of the background assuming infinite statistics 

The COSY-11 facility used to investigate the quasi-free pn pnX reactions 
is equipped with detectors for the registration of fast charged particles [18, 19], 
neutrons [20,21], and slow spectator protons [9]. Experiments are conducted using the 
COSY proton beam [22,23] incident on the internal deuteron cluster target [24,25]. 

The total CM. energy (y^ ) available for the proton-neutron scattering is fixed 
by the four- vectors of the incident and spectator protons. Due to the internal motion 
of the nucleons in the deuteron, this varies from event to event and in the rj region a 
range of about 60 MeV is scanned for a fixed beam momentum [14,26,27]. The four- 
momentum vectors of the beam and spectator protons are known for each registered 
event. Therefore, the collected data can be grouped according to the excess energy 
with respect to the pn pnr] process. The production of the r] meson can only occur 
if the excess energy Q, 

Q — y/s '^proton '^neutron ^77 (f) 

is positive. For negative Q only pions can be created. 

Fig. 1 shows in a schematic manner the mass distributions expected at three 
different values of Q for the unobserved system X produced via the pn — > pnX 
reaction. The upper panel presents an outlined missing mass spectrum when Q < 0. 
The shape of this continuous distribution arises from a convolution of the real invariant 
mass distribution of the multi-pion system, the efficiency of the detection setup for 
the simultaneous registration of nucleons from the pn pn pions reaction, and also 
from the experimental resolution of the missing mass determination. For Q > a 
signal from the r] meson is expected on top of the multi-pion mass distribution at a 
position corresponding to its mass, as shown in the middle panel of Fig. 1. For data at 
higher Q the peak from the rj meson will remain at the same position on the mass axis, 
though the distribution of multi-pion invariant mass will spread towards larger masses 
(lowest panel of Fig. 1), extending up to the kinematical limit. However, in practice, 
as will be argued later, the change of the multi-pion background can be approximated 
as a shift of the entire spectrum in the direction of higher masses. 

In the case of the pn pnrj reaction measured at the COSY-11 facility, the 
experimental mass resolution amounted to about 5 MeV(FWHM) [29]. Therefore, for 
measurements at only few MeV above the threshold the shape of the pion background 
must be well known in order to extract the 77 signal. The lack of a clear indication 
of the upper end of the signal resulting from the rj meson production (middle panel 
of Fig. 1) constitutes the main problem in disentangling the contributions from the r] 
meson and pions. 

For simplicity of argumentation, let us first assume in subtracting the multi-pion 
background that the number of registered events is very large, such that one can 
neglect any statistical fluctuations, and that the shape of the reconstructed invariant 
mass of the multi-pion background is independent of the excess energy Q. The form 
of the background could then be determined from a missing mass spectrum dN^ /dm 
from an infinitesimal range at any negative values Q. Having postulated that the 
shape of the background dN'^{m,Q)/dm does not depend on Q, it is convenient to 
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Figure 1. Schematic depiction of the missing mass spectra for the pn — » 
pnX reaction below {Q = Qi) and above {Q = Q2 and Q = Q3) the 
threshold for the rj meson production. The broadening of a signal from the 
7j meson with the increasing of Q is kinematical effect discussed in detail 
in reference [28]. 

express it in terms of a normalised function B{mrj + Q — m) of the difference between 
the kinematical limit (m^ + Q) and the given mass m. 

Being not limited by statistics, one could divide the positive Q range into very 
narrow subranges within which the resultant missing mass spectrum would be the 
simple sum of a signal from the rj meson and the form B multiplied by an appropriate 
factor. The discussed situation is depicted schematically in Fig. 2, where the lower 
panel shows the method of how to construct the background. 

If the two assumptions mentioned above were valid then, in order to derive a signal 
of the rj meson from a missing mass spectrum for positive Q, it would be sufficient 
to subtract from this spectrum a missing mass distribution determined for negative 
Q after the shift of the latter to the kinematical limit (dotted line) and normalisation 
at the very low mass values where no events from the 77 meson production arc 
expected (dash-dotted line). In such a case the contribution of the pn pnr] reaction 
could be extracted without the necessity of any additional assumption regarding the 
unknown distribution of the background expressed as a function of the excess energy 
dN'^/dQ. In this case the analysis would not demand an exact knowledge of the 
dynamics of the background reactions. 
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Figure 2. a) Distribution of the excess energy with respect to the 
pn pnri reaction plotted schematically for a beam momentum of 
Pi, = 2.075 GeV/c. The shape results from i) the genuine excitation 
functions for ij meson and multi-pion production in the proton-neutron 
collisions, ii) the distribution of the Fermi momentum of the nucleons in 
the deuterium target, and iii) the acceptance and efficiency of the detection 
system. For a detailed discussion of this issue see e.g. ref. [14,26]. The 
contribution to the spectrum from pion production dN^ /dQ cannot be 
derived unambiguously from this spectrum. Therefore, entirely arbitrary, 
two different possibilities are pointed out as dotted and dashed curves, 
b), c) Form of the missing mass spectra as derived for the negative {Qb) 
and positive {Qri) values of Q. To each value of Q a continuous spectrum 
of masses caii be assigned, originating e.g. from the pn — > pmrix reaction. 
However, there is always a unique relation between the value of Q and the 
maximum mass which can be created {rrimax = rrin + Q). Therefore, 
for better visualisation, the axis of Q for the dN/dQ distribution and the 
mass Eixis of the missing mass spectra dN/dm were arranged such that a 
maximum mass achievable for a given value of Q lies exactly below the 
corresponding Q value in the dN/dQ plot. 

d) The shape of the background in the mass spectrum for Q > 
can be constructed from the shape of the multi-pion mass distribution 
(shaded histogram) after its shift to the kinematical limit (dotted line) and 
subsequent normalisation at low masses (dash-dotted line). 
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3. Background estimation under the real experimental conditions 

Due to the finite statistics, one needs to study the missing mass spectrum in finite 
intervals of Q. This does not influence the shape of the missing mass distribution 
arising from the pn pnrj reaction, though it alters the form of the multi-pion 
background. The latter dN"" / dm{m) , when derived from the finite AQ range, becomes 
the convolution of the B function and the number of background events expressed as 
a function of the excess energy dN'^ /dQ. For the range spanned between Qmin and 
Qmax the missing mass distribution of the background can be expressed as: 



In the dN/dQ spectrum we cannot extract the function dN'^/dQ in a model- 
independent way for positive values of Q (compare dashed and dotted lines in Fig. 2a). 
We must therefore take as narrow ranges of AQ as reasonable, taking into account the 
experimental resolution in Q, and to assume that within such a range the dN'^/dQ 
function is constant. It is natural to take the width of the Q subranges equal to 
the FWHM of the resolution of the Q determination. Such a width choice partially 
validates the constant dN^dQ assumption, due to the unavoidable experimental 
smearing within the bin. 

By reviewing the experimental distribution dN/dQ shown in the upper panel of 
Fig. 3, one recognizes that the obtained statistics are indeed insufficient to derive the 
dN'^ / dm spectrum from the bin of Q with a width equal to the experimental resolution 
(i.e. FWHM = 5 MeV). The statistics could, however, be improved significantly if all 
events registered with negative Q could be taken into account. This can be realized by 
adding to the missing mass calculated for a given event an amount (Q^ — Qt), which 
will shift the background events measured at Qi, such as if they were determined at 
the excess energy of Qn- The resultant modified missing mass distribution obtained 
from the entire data sample of negative Q values can be identified with the function 
B{mri + Qn ~ ™) needed for the derivation of the background distribution within 
a finite excess energy bin centered around the Q,, value. In such a way, when 
investigating the missing mass spectrum for the Q,, bin, one would obtain directly 
the background spectrum which, after the normalisation, could be subtracted from 
the distribution containing the signal from the r/ meson. When the number of events 
is small, this approach enables one to reduce significantly the statistical fluctuations 
in the background spectrum. However, as was already stressed, the entire procedure 
could only be justifled if the shape of the background did not alter within the discussed 
range of the excess energy Q. Arguments in favour of this assumption will be presented 
in the subsequent section. 

4. Why does the experimental multi-pion background not change 
significantly with variation of the excess energy? 

We now present arguments supporting the supposition that the measured shape of the 
background does not alter significantly with changes in excess energy, at least over the 
energy range studied. It is important to stress that this assumption worked correctly 
when evaluating the data of the reactions pp — > pprj' [32] and pp pprj [28] and 
that the derived cross sections are consistent with those obtained at the SATURNE 





facility [33]. 
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Figure 3. Upper panel: Experimental distribution of the excess energy Q as 
measured with the COSY-11 setup, using a deuterium target and proton beam 
with a momentum of 2.075 GeV/c [14]. 

Lower panels: Scheme of the missing mass distribution for two arbitrarily chosen 
values of Q < 0. 



As already mentioned, the shape of the background spectrum reflects i) the real 
distribution of the invariant mass of the multi-pion systems, ii) the acceptance of the 
detection setup, and iii) the c;xperimental resolution. The acceptance of the detection 
system, its efficiency and resolution do not change significantly over the considered 
variation of excess energy. On the level of a few per cent, this is true even in the case 
of the detection setups characterized by a very limited geometrical acceptance, such 
as e.g. the COSY-11 facility [34]. The shape of the genuine background distribution 
is the only significant factor that is relevant to our discussion. In practice, from the 
high statistics measurements of the pp — > ppr] reaction, we know that the shape of the 
background could be sufficiently well described assuming that it was due to two-pion 
creation only [35-37]. Adding into the simulations three- and four-pion events was 
superfluous, and did not improve the agreement with the experimental distributions 
signiflcantly. 
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Figure 4. (Upper panel) Missing mass with respect to the pn system 
calculated for the pn — > prnviv process at the threshold for the pn — » pnr) 
reaction. Solid line shows uniform phase-space distribution, dotted line 
presents its modification by the final state interaction (FSI) between proton 
and neutron, and daKlicd line represents spectrum obtained assuming that 
the pious are produced via the pn — > AA — » pnnn reaction chain. Dashed- 
dotted line shows the effects of both the proton-neutron FSI and the 
resonance production together. 

(Middle panel) Comparison between the distributions calculated regarding 
resonance production and FSI for = 20 MeV (solid line) and 

Qr, = MeV (dotted line). The latter is shifted by 20 MeV. 
(Lower panel) Difference between the spectra of the middle panel 
normalised to the solid line: 

{dN''{Qr, = 20)/dm - diV'^(Q^ = 0)/dm) / {dN^'^Qr, = 20) /dm). 

Fig. 4 demonstrates the missing mass distribution of tttt system produced via 
the pn — > pnX reaction calculated taking into account experimental resolution of the 
detection system and various creation scenarios. 
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Missing mass [ GeV/c^ ] 




0,3 0,4 0,5 0.6 0,7 0,8 0,9 

Missing mass [ GeV/c^ ] 

Figure 5. (Upper panel) Missing mass with respcet to the pn, system 
calculated for the pn —^ pnnn process at the threshold (dotted Une) and 20 
MeV above the threshold (solid line) of pn — » pnr)' reaction. The dotted 
histogram was shifted by 20 MeV towards the larger masses. (Lower panel) 
Difference between the spectra of the upper panel normalised to the solid 
line: {dN''{Q^, = 20) /dm - dN^iQ,,, = 0)/dm) / {dN''{Qn> = 20) /dm). 

When Q increases the distribution of the multi-pion invariant mass broadens 
towards larger masses. In practice, if the change in the vahic of Q is small compared 
to the excess energy with respect to threshold of the pn pn pions reaction, the shape 
of the mass distribution at its edge remains approximately unaltered. Of interest for 
us is the range of about 40 MeV/c^ from the higher mass limit. The events represented 
by the solid and dotted lines in the middle panel of Fig. 4 were simulated at excess 
energies of Q = 20 and Q = MeV with respect to the pnr] system. Prom the 
comparison of the form of both spectra one realises that at the edges there is no 
noticeable difference between the shape of the curves. Thus, this shape, convoluted 
with an acceptance function independent of the Q value, will lead to the form of the 
background which will be a function only of the distance from the kinematical limit, 
over at least a range of about 40 MeV. 

The quantitative appraisals of the systematical error caused by the subtraction 
of the shifted and normalized spectra is depicted in the lower panel of figure 4. This 
picture demonstrates a fractional systematical error due to the discussed method. It 
was calculated as a difference between the mass distributions dN/ dm" at Q = 20 MeV 
and Q = MeV normalized to the dN/dm'"{Q = 20 MeV). From the figure one can 
infer that the systematical error in the missing mass range pertinent to our discussion 
(close to the kinematical limit) is in the order of a few per cent of the background 
value. In the middle and lower panels of Fig. 4 wc show the comparison for the 
production model depicted by the dashed-dotted line in the upper panel, however, we 
have checked that for other discussed possibilities the determined systematical error 
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is similar. 

It is worth noting that in the case of the 77' meson, due to its twice larger mass, 
the relative changes between dotted and solid hnes are even smaller (see upper panel 
of Fig. 5), and hence the advocated assumption is even better justified. The lower 
panel of Fig. 5 indicates that the fractional error due to the changes in the shape is 
in the order of one per cent only within the range of 200 MeV in the vicinity of the 
kinematical limit. 



5. Example of application 

Recently the 77 and 77' mesons production in the proton-neutron collisions has been 
studied close to the kinematical threshold at the COSY-11 facility [14,15,17]. Although 
the evaluation of these data is still in progress, the first trials to distinguish between 
the multi-pion background and the rj meson production in the missing mass spectra 
lead to a clear signal for the pn pnr] reaction. Fig. 6a shows the sum of the 
missing mass spectra for all Q > bins (solid line) and the sum of the corresponding 
background spectra (shaded histogram). 




O 

O 20 



missing mass [ GeV / c ] 



a) 



0.35 0.4 0.45 0.5 0.55 
missing mass [GeV/c^] 



b) 



Figure 6. a) The solid line shows the missing mass distribution of the pn — > pnX 
process determined for Q > with respect to the pn — > pnrj reaction. The sum 
for all AQ intervals is shown. The shaded histogram depicts the missing mass 
spectrum determined for Q < 0. The procedure of the background construction 
is as described in the text, b) The points denote the experimental missing 
mass spectrum for Q > after subtraction of the multi-pion background. The 
superimposed solid line, normalised in amplitude to the data points, results from 
a Monte-Carlo simulation. 



The data were evaluated as follows. First the sample of events with excess energy 
greater than zero was divided into a few intervals of AQ — 5 MeV and a missing 
mass spectrum was calculated for each sub-sample separately. Next, for each of 
these spectra, determined in interval AQ around a given Q,, value, a background 
was constructed from all events with Q less than zero. The missing mass of each 
background event, measured at given excess energy Qfci was increased by an amount 
{Qri — Qb), as depicted in Fig. 3 and described in section 3. For each AQ interval 
separately, the background spectrum was normalised to the distribution at a given 
for mass values less than 0.25 GeV/c^. In this region the events correspond to a single 
pion production, and since the energy is high above the threshold the cross section for 
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TT^ production stays nearly constant when the excess energy changes only by few tens 
of MeV. 

In order to make a check of the background extraction we have constructed 
it dividing the range of Q values into two parts equal in number of counts. The 

comparison of the spectra resulted from the two parts is shown in Fig. 7. Indeed, 
the determined spectra from different ranges of Q are consistent within the statistical 
errors. 

I 100 

3 
O 

o 

60 



20 

0.4 0.44 0.48 0.52 0.56 

Missing mass [GeV/c^] 

Figure 7. Background distribution obtained using events with Q ranging 
from -9 MeV to MeV (dotted line) and from -40 MeV to -9 MeV (solid 
line) . 

Fig. 6b depicts the spectrum of the missing mass for Q > after the background 
subtraction. The superimposed line indicates the expectations from the Monte- 
Carlo simulations, taking into accomit the momentum spread of the beam [39], 
Fermi momentum of the nucleons inside the deuteron [29-31], beam and target 
dimensions [39], multiple scattering in the detector components and air, and position-, 
time- and energy- resolution of all detector components [18-20,35,40]. Clearly, within 
the statistical uncertainties, the results of simulations conform very well to the shape 

of the observed distribution. 

For a better visualisation, the sums for all AQ intervals arc shown in Figs. 6a 
and 6b because the individual spectra arc statistically much less significant, due to the 
small number of events in this first measurement, which was meant only as a feasibility 
test of the COSY- 11 facility for the study of the meson production via quasi- free 
proton-neutron scattering. A more detailed view of the region around the r] peak is 
shown in Fig. 8a. For comparison we present also in Fig. 8b the result obtained before 
the current method was elaborated [14]. This first analysis was done by deriving a 
missing mass background spectrum from all events with Q < without adding (event- 
by-event) the difference of excess energies between a corresponding background event 
and that from the bin for which the signal of r] meson was being studied. 
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Figure 8. Missing mass spectra: a) determined according to the method 
described in this article, and b) using an earlier approximate procedure [14] . The 
solid line depicts results of the simulation normalised to the data. 



Although at first sight both spectra look similar, on closer inspection one 
recognizes that the left hand histogram agrees statistically much better with the 
experimentally extracted data than that on the right. The difference is understandable 
since the shift of the background to the kinematical limit diminishes the number of 
background events on the left of the eta-peak and increases that on the right (see 
Fig. 9), leading to the better agreement with the shape predicted by the Monte- 
Carlo studies. Moreover, the total number of events assigned to the pn — > pnrj 
reaction, when extracted according to the new procedure, is by about 30% lower 
than previously estimated [17] and now conforms better to the results obtained by 
the WASA/PROMICE collaboration [6], where photons from the decay of the 77 were 
detected. 
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Figure 9. Form of the background determined according to the method 
described in this article (dotted line), and using an earlier approximate 
procedure [14] (solid line). 

Finally in Fig. 10 we present a dN^ /dQ function determined for the positive 
values of Q. One can see that the distribution dN'^ /dQ for positive values of Q is rather 
not constant within the AQ bin contrary to the assumption made when constructing 
the background. Therefore, we have made an iterative procedure and construct a 
background taking in formula 2 the obtained distributuion of dN^ /dQ. The result 
was statistically the same even after the first iteration. This can be understood since 
the width of the AQ bin is rather small compared to the mass resolution. One bin at 
the mass spectrum (eg. in figure 8 ) corresponds to 4 MeV. Therefore, even though 
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the dN'^/dQ significantly differs from being constant it hardly influences the dN/dm 
spectrum within the present statistics. A larger sample of data will enable to determine 
dN'^/dQ iteratively more precisely and hence will permit to obtain results not biased 
by the assumption that dN'^/dQ is constant within a studied bin. 



Experiment 
Extracted background 




20 40 

Q [MeV ] 



Figure 10. The solid line shows experimental distribution of the excess 
energy Q as measured with the COSY-11 detection setup (see also 
figure 3), and the superimposed dotted histogram indicates a spectrum 
of (IN^ / dQ(Q) extracted for positive values of Q taking into account all 
events of Q < and the method introduced in this article. 

6. Summary 



A technique for disentangling contributions of multi-pion and single meson production 
in the missing mass spectrum of quasi-free pn pnX reactions has been presented. 
The method enables one to extract a signal for the pn pn meson reactions measured 
close to the threshold, when the mass resolution is comparable with the centre-of-mass 
excess energy. In this energy range the peak in the missing mass spectrum from the 
single meson creation is close to the kinematical limit and the separation of the signal 
due to the single meson from that of the continuous distribution associated with 
multi-pion production becomes nontrivial. It was also shown how to combine events 
of different excess energies in order to gain statistical significance in the determination 
of the shape of the multi-pion mass distribution. 

In the case of the pd — > pnXp^p reaction measured by the COSY-11 
collaboration [14], we demonstrated that the method leads to the derivation of the 
signal from the quasi-free pn — > pmr} process that is consistent with expectations. This 
successful application, and the fact that the technique is independent of the detection 
system, suggests that it can be applied more generally for distinguishing between 
single and multi-meson production in missing mass spectra obtained in pd pnXpsp 
reactions. 

The presented method will be used for the extraction of the signal from data on 
T]' meson production via the pn — > pnr]' reaction [41], whose rate is expected to be 
about forty times lower than that for the r] meson. These cross sections are utterly 
unknown and their determination may be of crucial importance for the study of the 
gluonic degrees of freedom in the 77 and rj' meson production processes [42]. 
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